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Abstract 

We study dilute gases with short range interactions and large two-body scattering lengths. At 
temperatures between the condensation temperature and the scale set by the range of the potential 
there is a high degree of universality. The first two terms in the expansion of thermodynamic 
functions in powers of the fugacity z, which measures the diluteness of the system, are determined 
by the scattering length only. The term proportional to z 3 depends only on one new parameter 
describing the three-body physics. We compute the third term of the expansion and show that, 
for many values of this new parameter, the z 3 term may be the dominant one. 
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I. INTRODUCTION 



Only on special circumstances the thermodynamic functions of a system can be evalu- 
ated starting from the microscopic interactions. However, in many situations only a few 
characteristics of the microscopic interactions are relevant. All systems sharing these same 
microscopic characteristics have then the same macroscopic behavior, that is, there is a cer- 
tain degree of universality. An example of such a system is a gas of particles with short range 
interactions. As long as the density and temperature are such that the typical wavelength 
A of the particles is much larger than the range R of the forces, the details of the interac- 
tion potential is largely irrelevant and the thermodynamics is determined by the two-body 
scattering length a, up to corrections of order X/R. This observation has been used both for 
fermionic and bosonic gases since the 1950's |3|, f|, |^, [], 0, [5| in the computation of 

dilute gas properties in an expansion in powers of na 3 , where n is the density of particles. 

Typically, the size of the two-body scattering length is comparable to the range of the 
force, a ~ R. However, there are important situations where the interactions are fine-tuned 
in such a way as to make a 3> R, even though other low-energy parameters like the effective 
range ro, etc., still have the size expected on dimensional grounds ro ~ R. We have in 
mind two of these situations. The first one, a gas of neutral atoms, has received enormous 
experimental and theoretical attention recently. The range of the interactions between the 
atoms is set by the length scale R ~ (MC^/fi 2 ) 1 ^ of the van der Wals force — Ce/r 6 . In some 
atomic species like 87 Rb and 4 He the fine tuning needed for large values of a is provided by 
Nature. In the case of 4 He atoms for instance, a = 104v4, a value much larger than R ~ 5A 
or the effective range r ~ 7 A. More importantly, an external magnetic field can be applied 
in order to artificially modify the scattering length, making it a tunable parameter (Feshbach 
resonance) . The second example of unnaturally large scattering lengths is a gas of neutrons 
(and protons, if the Coulomb force can be disregarded). The range of the nuclear forces is of 
the order of the Compton wavelength of the pion (R ~ h/m^ ~ 1.5/m), but the scattering 
lengths are significantly larger (5.42 fm for the proton-neutron in the spin-triplet state to 
—18.8 fm for neutron-neutron in the spin-singlet state). 

The typical momentum of the particles in a gas is set by the larger of the inverse inter- 
particle distance n 1//3 and the inverse thermal wavelength A = (MT/2tt)~ 1 ^ 2 (from now 
on we will use h = 1). In the natural CcLSG, CI « R, the universal regime occurs for 
1/A, n 1 / 3 < 1/a m 1/R. This regime is essentially perturbative in the sense that, at any 
given order of the expansion in powers of 1/R, a properly set up diagrammatic expansion 
involves only a finite number of diagrams jlOfl . On the other hand when n 1 ^ 3 « 1/R 
the details of the particle interactions are important and each system should be studied in 
a case by case basis. The presence of large scattering lengths opens up an intermediate 
regime, where the typical wavelength A of the particles is comparable to a, but still smaller 
than R . The problem is no longer perturbative, as evidenced by the fact that bound states 
of size ~ a are expected, but some degree of universality should still hold. This regime 
can be attained at very low temperatures and a moderately high densities na 3 ~ 1 ^ nR 3 
or at a moderately high temperatures a/ X ~ I > R/X and a very low densities. It is an 
outstanding problem to understand the first of these regimes, since, most likely, many-body 
correlations are not suppressed, and a number of publications have appeared recently on the 



subject [IT, 0, [13], |I3, |I^, In this paper we will consider the second case, namely a 
dilute (na 3 < 1), moderately hot ((MT ' /2-n) 1 ^ 2 a fa 1) gas of resonating a ^ R particles. 
Since we are considering low densities, it is convenient to phrase our discussion in terms 
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of the virial expansion, where the different thermodynamic functions are expressed as power 
series on the fugacity z = e^, where T = 1/(3 is the temperature and \x the chemical 
potential. The particle density n and the pressure P are given by 

n = + 2b 2 z 2 + 3b 3 z 3 + ■■■), 

A 

P = ^(b lZ + b 2 z 2 + b 3 z 3 + ---) . (1) 

The usefulness of the virial expansion resides on the fact that, for small densities, n < A -3 , 
z is also small z <C 1. The coefficients bi contains contributions coming from m-particle 
correlations for all m < I so its computation involves the solution of the /-body problem. 
The calculations of b\ and b 2 are rather easy and can be done analytically in the system 
considered here. The computation of b 3 is a little more involved: it includes some numerical 
integrations, and it is related to the quite unusual properties of the three resonating particle 
system. 



There is a general formula relating 63 to S-matrix elements [fTfl| . However, the relevant 
matrix elements relevant here are the ones describing a variety of processes like three-to- 
three-particle scattering, dimer break up in a collision to a particle, etc.. It seems to us 
that the equivalent method followed below implicitly includes all these processes in a simple 
way. Also, since it is an straightforward application of standard effective field theory and 
many-body physics methods, it may have some methodological interest. 

We will use the language of effective field theory, which is natural when exploring univer- 
sal, low energy, large distance properties that are independent of the short distance details. 
We will work at the leading order in the small momentum expansion where the Lagrangian 
of the system is 

where ip (ip^ 1 ) is the field that annihilates (creates) a particle, Co and Do are constants that 
will be determined later and the dots stand for terms with either more derivatives or ip 
fields, whose contributions are suppressed at the order we are considering here. We will 
discuss the bosonic case first and comment on the few, but important, differences present 
in the fermionic case later. It is convenient to introduce a dummy field d with the quantum 
numbers of two particles (a dimer) and use the equivalent Lagrangian 



C = ^{id + ^)V> + Arfd - ^(^V 2 + h.c.) - g 3 dU^iP + • • ■ , (3) 

where 2g 2 /A = Co and 36g 3 g 2 /A 2 = D . The equivalence between the two Lagrangian 
can be seen by performing the gaussian integration over the auxiliary field d and recovering 
Eq. (H). The value of the constant A is arbitrary and affects only the normalization of the 
dimer field. Physical quantities depend on it only through the combinations Co and Do- 

We now compute the particle density n in a expansion in powers of z and, by comparing 
with Eq. ([[]), determine the coefficients b\. The computation of the two first virial coefficients 
is rather trivial, and there are fairly explicit general formulae for them. We will quickly 
discuss them here in order to explain the method used in selectingwchich diagrams contribute 
at each order. 



3 



FIG. 1: Graph determining the density at leading order in z. The cross represents an insertion of 
the particle number operator. 



In selecting the diagrams contributing to given order in z one should notice that diagrams 
with a closed particle line vanish in vacuum (z — 0). Consequently, at low densities their 
contribution is suppressed by one power of z for each closed particle loop. With that 
observation in mind we see that the only diagram contributing to the particle density at 
leading order in z is the loop diagram shown in Fig. ([!]): 



k 



d 3 k -1 



(2tt) 3 -ik + e k 

z 2 z 3 



The sum over the frequencies ko is over all integer multiples of 2ttT and e k = k 2 /2M — fi. 
Eq. (|4]) determines b\ — 1, which is the free gas result. It also gives some contributions to 
6 2 and 6 3 . 

The 0(z 2 ) contributions can be divided into the one-body contribution coming from 
the second term in Eq. @ and the two-body contributions . From Eq. (||) and Eq. ([[]) 

(1) 5 (2) /n=3 

we find V£' = 2-1. The dia grams contributing to n 2 are shown in Fig. (y). The need for 
the resummation in the full dimer propagator indicated at the bottom of Fig. (0) is more 
easily explained after computing it. The dimer propagator is given by a geometrical sum 

p(p) = -l + i EW i + - 

(5) 



A + E(p)' 
where 

E(p) = ^\p) + rt l \p) + 0(z 2 ) (6) 



Mg 2 [2k p 2 ^ \ zMg 2 [°° 7 2 0^ pk 

—r-\ \ ~. — 2Mu-tMp ir- / dfcV^arctghf 

4tt I 7T V4 >* 2vr 2 p 7 * y k 2 + ^-2M^-iMp Q 



In particular, the leading order dimer propagator is the same as in the vacuum 

4tt 1 
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FIG. 2: Diagram determining the density at leading order in z 2 (top) and the dressing of the 
dimer propagator. 



The particle-particle scattering amplitude in the vacuum is given by g 2 T>. Thus, by setting 
AnA/Mg 2 + 2 A/ 7i = l/a, we obtain the leading order scattering amplitude in the effective 
range expansion (after analytically continuing the amplitude to real values of the energy by 
setting ko — > k — 2M/J, + it) . Notice the fine tuning between the linearly divergent term 
2A/V (kinetic energy) and the interaction term 47rA/Mg(A) 2 (potential energy) needed to 
produce a small value of l/a. We can understand now the reason for the resummation of 
graphs at the bottom of Fig. (Q). For momenta of order p ~ y/MT ~ l/a, the square root in 
the denominator of Eq. (0), generated by the loops in the dimer propagator, is not negligible 
compared to l/a. 

Plugging the expression in Eq. (|7|) in the diagram in Fig. (0) we have 



MTJ2 J 



d 3 k 1 1 



k J ( 27r )"-I + ^ -2Mfi-iMk yJ^ -2M/i-iMk 

, r . dri f d 3 k 1 1 1 

MT A / , 7= • (8) 



2vu J (2vr) 3 e^ - 1 i j& 



(2) 



The integration over 7], on a contour encircling the imaginary axis, can be deformed into an 
integral over the discontinuity on the branch cut (which describes the dimer break up) plus 
the contribution of the dimer pole. The final result is 

^^^^^(uMy), ( 9) 

where B2 = l/(Ma 2 ) is the location of the pole in the dimer propagator. This pole corre- 
sponds to a bound state (if a is positive) or a "virtual" bound state, that is, a pole in the 
unphysical sheet (if a is negative). Eq. (B|) is a particular case of the classic formula relating 
b 2 to the two-particle scattering phase shift ]T8| 



6 <W/ 2(e ^ + i p^ e ^) , (10) 



in the case where the phase shifts are given by the leading order effective range expansion 
k cot 5(k) = — l/a. A few points are worth mentioning here. First, 62 is a continuous 
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FIG. 3: Graph determining the density at leading order in z 3 and the integral equation for the 
three-body amplitude. 

function of a at 1/a = 0. The change in the potential needed to take a from a large and 
positive value to a large and negative value is small, and that change has only a small 
effect in the thermodynamics, even though there is a qualitative difference in the spectrum 
(from a bound state to a virtual bound state). Second, b 2 is positive and, for B 2 > T, it 
is exponentially large. In the case a > 0, this is easily understood as a consequence of the 
existence of a bound state: when TcB 2 most particles form tight two-body bound states 
therefore increasing the density for a fixed temperature and chemical potential. In fact, if 
the z 2 terms dominate, the ratio of the two equations in Eq. (P give the equation of state 
of a free gas with density n/2. 

The third virial coefficient 63 includes three pieces. The first, 63 , is determined by the 



piece in Eq. (HI) and equals b 



(i) 



The second one, b^\ comes from a subleading 



piece of the diagram on Fig. (En and is given by 



n 



(2) 



3 b 



(2) 



(11) 



2T£ 

fco 
fco 

i(3z 3 r°° 



d 3 k 



(2vr)3 

d?k d 
(2tt) 3 ^ 



(z><°> (/c)) 2 s (i) (k) 4i-x (o) (k) + p (o) (k) o) 



(P (0) (A;)S (1) (A;)) 



7T 4 , 

7T ./ 



<//. " / dp I dQe /3 2Af a rctgh( 



p k 



k 2 + p T - 2M/i - iM Po _l + Jk _ MC 



where the integral over ( runs on a vertical line just to the left of the pole at — B 2 . The 

integral in Eq. (|TT|) cannot be expressed analytically and has to be computed numerically. 

(3) 

The third piece, b 3 , comes from actual three-body correlations in the system, and brings 
in new physics besides that contained in two-body scattering. The diagrams contributing 
to it are shown on Fig. (fj). Just as the two-body scattering amplitude cannot be computed 
in perturbation theory, the three-body amplitude that enters in Fig. (M also involves a 
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resummation of an infinite number of diagrams. This is because each additional loop involves 
a factor of pa, where p is the typical momentum flowing through the loop. In the case 
considered here Q~l/A~l/a and the "suppression" factor is of order 1. The diagrams in 
Fig. (§) gives 



"P = T ^JW?W? (^) ><P + *)n-V,V + k;- P , P + k) (12) 

= E S§?W?^^ (^) 2 + + * -p. p + *) 

where in the second line p = (— z?7, p) and T(—p, p+k; —p, p+k) is the forward particle-dimer 
scattering amplitude determined by the diagrams at the bottom of Fig. (^). The integral 
over T) is dominated by the particle pole at r\ = —e p , with the contribution coming from the 
dimer pole and cut suppressed by two powers of z. We then have 



2ttz 3 f d 3 k d 3 P f d( ar rrl „ ifc 2k n k 2Jfe, 



= , (13) 



where we use the center-of-mass variable P = (ie p ,p + k/3). The integration over the new 
variable ( = 3/i + ( in on a vertical line just to the left of the left-most singularity that can 
be a trimer pole at k 2 /QM — B% or a dimer pole at k 2 /6M + 3P 2 /4M — I^- The cuts in the 
three-body amplitude describing dimer-particle scattering, dimer breakup, etc., are on the 
real axis, to the right of 2 and 3 particle poles. 

Unlike the dimer propagator, the diagrams adding up to the three-body amplitude 
T(— p, p + k; —p,p + k) do not form a simple geometrical series and cannot be summed 
up analytically. However, their sum is determined by the integral equation depicted at the 
bottom of Fig. (|j). This integral equation (Faddeev equation) is particularly simple for the 
separable potential used here 1 and was derived many times before with a variety of tech- 



niques p0| , 21, ^2|. It is more easily written in terms of the s— wave amplitude P, P') 
defined by 

Mg 2 g %(C,P,n = f dP dP' T(-P + §, P + f ; -P + §, P + f ) | fc0= ^-< > (M) 



P 2_4M (c + j g 2) J 4n 4n _i, A /^T,£ 



Higher partial waves, describing more peripheral particle-dimer collisions, give suppressed 
contributions. %((, P, P') satisfies 

T (C, P, P') = K(C, P, P') + - T dqq 2 ^%^—K(C, q, P'), (15) 

n Jo <r - -»-(C + B 2 ) 



1 Some simplifications in the computation of 63 in the case of separable potentials were noted in |p_9|| 
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with the kernel 



K((,P,P') 




PP 



7 log 



P 2 + pp, + p,2 _ M £\ g 3 ( A ) 



p2 _ pp, + p/2 _ M £ 



M# 2 



(16) 



Eq. ( |T5D has a number of surprising properties [f23], |24], |25|, [2"T |. 
more naturally described in terms of renormalization theory. It was shown in 
To(C,P,P') is kept cutoff independent for small values of P,P' < A if, and only if, the 
three-body force varies with the cutoff A as 



This properties are 
@ £3 that 



#3(A)A 2 = sinQo log(A/A*) - arctan s ) + ^ . J_, 
2Mg 2 sin(so log(A/A*) + arctan Sq) aA' 



(17) 



where A* is a new parameter that cannot be measured through two-particle experiments 
and So ~ 1.006 is the solution of a certain trigonometric equation. The failure to include the 
three-body force term leads to a strong cutoff dependence on the results and thus, to the 
impossibility of arriving at model independent results. The three-body scattering amplitude 
depends, even at the leading order in the expansion in powers of X/R, on the value of the 
three-body force, here parameterized by A*. The parameter A* encodes all short distance 
physics, besides the value of a, necessary to describe three-body systems. Like a, it varies 
from one system to another and only the measurement of some three-body observable allows 
us to fix it. The appearance of an extra parameter not determined by two-body scattering 
precludes the possibility of predictions in the three-body sector based only on the value of 
a. 

Using Eq. (|l^) and Eq. ([14]) and performing a trivial integral over k and the angular part 
of P we arrive at 



n 



(3) 



3^ 

3 ^3 



.V3 



4vr 2 A d 



/ dPP 2 d(e 



%(C,P,P) 



' 3 -f-M(-B 2 



l + i 
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The amplitude T ((, P, P) and the quadratures in Eq. (|) have to be performed numerically. 
This is not computationally demanding and a Mathematica notebook that performs this 
task can be downloaded from [http : //www-nsdth . lbl . gov/bedaque| . 



The coefficient 63 = + + bf^ is a function of the temperature and of the pa- 



rameters describing the microscopic dynamics (the scattering length a and the three-body 
parameter A*). In Fig. ([D we show the value of 63 as a function of 1/a for a few values 
of these parameters. As an example we pick M = M4 He , T = 10~ 7 eV ~ 1.16mK and 
-A 2 g 3 (A)/(2Mg 2 ) = H(A) = -3.235 (lower blue solid curve), -0.22 (intermediate red solid 
curve) and 0.05 (upper black solid curve). In all of them we use A = 200 eV, but, as men- 
tioned above, the results are cutoff independent up to small corrections of order 1/Aa. The 
first of these values was chosen following [JJIJ and it gives rise to a trimer with the bind- 
ing energy of the shallower 4 He trimer (as predicted through potential model calculations 
[31], |33|] ) when a is set to the value inferred from the dimer measurement a = lO'tfgA 
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For the value of the A used the deeper trimer state is absent. 
When a and A* are such that there is a three-body bound state B 3 > B 2 , we find 
empirically that 63 scales roughly as e^ B3 , see Fig. |j. This is not surprising and is a direct 
analogue of the physics described by Eq. (P) in the case of the coefficient b 2 : when j3B 3 3> 
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b 3 



1000 




1 



20 a 



FIG. 4: &3 as a function of the inverse scattering length for three values of H(K = 200eV) (solid 
lines, H = —3.23,-0.22,0.05, from bottom to top). The dashed line is the estimate e^ Bs . 



1 most particles form three-body bound states, which increases the density for a given 
temperature and chemical potential as compared to the free gas. This does not mean, 
however, that the integral in Eq. (|) is dominated by the trimer pole: the contribution from 
the two-and three-particle cuts is never negligible. For values of ah* for which B% approaches 
B 2 the dimer-particle cross section diverges. Still, the three-particle correlations described 
by 63 are smooth at those points. For systems with a trimer deeper than the dimer, 63 can 
be much larger that both b\ and 62, and dominate the virial expansion. This is not a rare 
situation, due to a well known but strange feature of the three-resonating particle system: 
as the potential is changed to make the dimer shallower , the trimer becomes deeper |E3 . 
This is what happens in the H = 0.05 case shown in the upper (black) solid curve in Fig. 

In atomic traps close to a Feshbach resonance it is a that is a tunable parameter. The 
influence of the magnetic field on the value of the three-body force parameter A* is small, 
since that describes short- distance physics of an energy scale much higher than the magnetic 
field, and can be disregarded in a first approximation. Thus we can regard the horizontal 
axis in Fig. (ffl) as denoting the magnetic field using the Feshbach resonance formula 



a(B) = a(B = 0) 1 



5B 



B — Bo 



(18) 



Let us summarize the approximation performed here that allowed such a simple evaluation 
of such a complex three-body dynamics. First, we approximate the two-body potential 
by a short range, delta-function-like interaction. Corrections to this approximation are 
proportional to the effective range r /X, where r ~ R is the effective range, much smaller 
than a by assumption. Effective range corrections can be easily included, even to very high 
orders, without spoiling the simplicity of the method, as it has been done in few-nucleon 
physics |2^. |35|] . The second approximation was the non-inclusion of higher partial waves. In 
the two-body sector, p-waves and higher are suppressed by at lest two powers of Rj X, and are 
small. In the three-body sector the suppression of the particle-dimer /-wave interactions are 
not suppressed by powers of R/X, but by 1/(1 + 1). All the higher partial waves correspond 
to a repulsive kernel in Eq. (|H>D, which do not support bound states and cannot produce an 
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enhancement of the form e 1333 . In practice the phase shifts are rather small but, if needed, 
they can also be easily included by solving the analogue of Eq. (|15|) corresponding to the 
higher partial wave and adding it to Eq. (pT3|). 

Our results are valid only in true thermal equilibrium, after two- and three-particle bound 
states had the time to form, and assuming that they stay in the system ( do not escape from 
the trap, if that is the case). For systems like the alkali atoms studied in magnetic/optical 
traps, that have deep bound states with typical interparticle distance of order ~ R that lie 
outside the validity range of our effective theory, this means that our results are relevant 
only for the metastable state before the collapse of the system, but after the two- and three- 
body bound states form. The formation of a /—body bound state requires the approach of 
/ + 1 particles in order to conserve energy and momentum, and their rates are consequently 
suppressed by n l+1 . These rates are not known at finite temperature but have been studied at 



zero temperature in [3]], |38| (recombination into shallow states) and |39| (recombination 
into deep states). In cases where a > the recombination rate into deep (two-body) bound 
states is estimated to be much smaller than the rate into shallow bound states suggesting 
that there is a time window in which our results apply. In the a < case there is no two- 
body bound state and it is not known how the rates for the formation of deep and shallow 
three-body bound states compare. 

Finally, let us consider the changes introduced in case of fermionic particles. In the non- 
degenerate case considered here, the effect of the quantum statistics in the thermodynamics 
is minor, and amounts to a change in sign in some of the coefficients b^' ^ ■ The elementary 
collisions, however, may differ a great deal due to the exclusion principle. If we have only 
one fermionic species in the system, s-wave scattering is impossible and all virial coefficients 
&2)&3, ■•■ are suppressed. In the case of two fermionic species (as a dilute neutron gas), 
two-body collisions are possible and the standard result in Eq. @ is valid. The physics of 
the three-body correlations is however, very different. No three-body force term without 
derivatives exist, and the three-body force contribution is suppressed. 63 can be computed 
in terms of a alone, but cannot ever be large and dominate the expansion, as the kernel 
appearing in Eq. (|T5|) would be repulsive and would not support a bound state. The case 
with three or more fermionic species with all the scattering lengths large but not necessarily 
equal, which includes the dilute nuclear matter case (protons and neutrons with spin either 
up or down), is very similar to the bosonic case. Three particles can occupy the same point 
in space without violating the exclusion principle and, as a consequence, the two coupled 
equations that substitute Eq. (|i5|) have very similar properties to the bosonic equation pOl . 
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